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SUMMARY 


This  report  concerns  the  problem  of  constructing  solutions  for 
transonic  flows  over  symmetric  airfoils.  The  aspect  of  the  problem 
emphasized  is,  of  necessity,  not  how  to  form  a  solution  for  compressible 
flow  but  how  to  simplify  the  initial  phase  of  the  problem,  namely,  the 
mapping  of  the  incompressible  flow.  In  the  case  of  the  symmetric 
Joukowski  airfoil  without  circulation,  the  mapping  is  relatively  simple, 
but  the  coefficients  in  the  power  series  are  difficult  to  evaluate.  As 
a  result,  the  problem  requires  simplification.  Instead  of  the  exact 
incompressible  flow  past  the  airfoil,  an  approximate  flow  is  used,  which 
is  derived  from  a  combination  of  source  and  sink.  This  flow  differs 
only  slightly  from  the  exact  one  when  the  thickness  is  small.  By  the 
same  method,  the  flow  with  circulation  is  also  considered. 

After  the  incompressible -flow  functions  are  approximated  in  this 
fashion,  the  numerical  calculation  of  the  corresponding  compressible  14 
flow,  by  the  hodograph  theory,  does  not  present  any  essential  difficulty. 


INTRODUCTION 


The  problem  of  transonic  potential  flows  has  been  considered  in  two 
previous  reports  (references  1  and  2).  The  object  was,  in  the  first 
place,  to  construct  a  solution  for  a  closed  body,  and,  secondly,  to  devise 
a  method  by  which  the  flow  can  be  easily  calculated.  The  method/  of  con¬ 
structing  a  solution  is  essentially  a  method  of  analytic  continuation. 

That  is,  a  complex  potential,  which  is  known  as  a  function  of  the  complex 
velocity,  is  represented  by  a  Taylor  series  about  the  origin  of  the  hodo¬ 
graph  plane  and  by  another  series  in  an  annular  region.  Both  series  are 
convergent  in  their  respective  domains  and  agree  on  the  common  circle  of 
convergence.  In  the  interior  region,  a  solution  for  compressible  fluid 
can  immediately  be  obtained  by  replacing  the  proper  particular  integrals 
in  the  Taylor  series.  The  new  solution  does  not  affect  the  radius  of 
convergence  and,  as  Mach  number  tends  to  zero,  reduces  to  that  of  the 
incompressible  flow.  The  solution  for  the  annular  region  cannot  be 
formed  by  simply  replacing  the  corresponding  particular  integrals  in 
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the  outside  series  of  the  incompressible  flow,  because  the  solution  so 
obtained  does  not  agree  with  the  inside  series  on  the  circle  of  con¬ 
vergence  and,  consequently,  does  not  represent  the  same  function. 

The  problem  is  then  not  so  simple  as  it  first  appears  to  be.  In 
order  to  make  both  series  represent  one  solution,  the  idea  of  retaining 
the  coefficients  in  both  series  of  the  incompressible  flow  must  be 
abandoned.  Instead,  both  the  coefficients  as  well  as  the  form  of  the 
solution  are  modified  to  satisfy  the  conditions  on  the  circle  of  con¬ 
vergence  (reference  2).  A  different  approach  to  the  same  problem  is 
that  due  to  Lighthill  and  Cherry  (references  3  and  4).  Instead  of 
choosing  in  advance  the  form  of  the  outside  series,  they  first  transform 
the  Taylor  series  into  a  double  series,  which,  on  interchanging  the 
order  of  summation,  can  be  summed  to  yield  a  single  series.  As  a  result, 
the  series  is  valid  anywhere  except  at  the  singularities  in  the  domain 
considered,  and  analytic  continuation  is  automatically  accomplished. 

The  continued  outside  series  consists  of  two  series;  one  is  derived 
from  the  outside  series  of  the  incompressible  flow  and  the  other  is  a 
Taylor  series. 

The  status  of  the  holograph  method  is  then  this:  If  a  complex 
potential  is  given  in  series  form,  a  corresponding  solution  for  com¬ 
pressible  fluid  can  be  constructed,  the  body  shape  being  unknown  until 
the  end  of  the  calculation.  It  is  apparent  that  the  fundamental  problem 
of  determining  the  flow  pattern  about  a  given  body  remains  unsolved. 
Moreover,  even  the  transformation  of  the  complex  potential  often  poses 
many  practical  difficulties.  Nevertheless,  the  holograph  method  is 
still  the  only  one  available  by  which  an  exact  solution  may  be  sought. 

For  this  reason,  there  is  still  ground  for  further  exploration  along 
this  line. 

It  is  the  object  of  this  report  to  examine  a  case  of  more  practical 
interest  than  that  previously  considered,  in  order  to  discover  possible 
means  of  simplification.  As  a  first  attempt,  a  symmetric  Joukowski  air¬ 
foil  is  studied.  In  this  case,  if  the  rule  of  transforming  incompressible 
to  compressible  flow  is  strictly  followed,  the  calculation  becomes 
extremely  laborious.  However,  if  a  special  case  of  small  thickness  - 
which  is  also  of  practical  interest  -  is-  considered,  the  functions  can 
be  simplified  to  such  an  extent  that  numerical  calculation  is  practical. 
The  process  is,  in  a  sense,  approximate.  But  as  far  as  the  compressible 
flow  is  concerned,  the  merit  of  a  method  lies  only  in  whether  it  yields 
a  good  aerodynamic  body  or  not .  It  has  never  been  demonstrated  that  an 
exact  incompressible  potential  function  gives  a  better  result.  From  this 
point  of  view,  this  is  just  another  way  of  choosing  a  complex  potential. 
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Finally,  it  must  be  emphasized  that  the  present  investigation  is 
only  the  beginning,  as  far  as  practical  calculations  for  useful  airfoil 
profiles  are  concerned.  Taking  the  approximate  function  as  a  guide, 
refinement  is  not  difficult  to  make  by  introducing  additional  terms 
according  to  the  character  of  the  function.  For  instance,  when  circu¬ 
lation  is  present,  the  hydrodynamic  functions  for  compressible  flow 
transformed  from  those  of  incompressible  flow  for  a  symmetric  body  will 
give  a  negatively  cambered  airfoil.  The  necessity  of  applying  correc¬ 
tion  terms  is  obvious. 

This  work  was  done  at  the  Graduate  School  of  Aeronautical  Engineering 
of  Cornell  University  under  the  sponsorship  and  with  the  financial  assist¬ 
ance  of  the  National  Advisory  Committee  for  Aeronautics. 


FUNDAMENTAL  EQUATIONS  AND  THEIR  PARTICULAR  SOLUTIONS 


A  steady,  potential,  and  isotropic  motion  of  a  perfect  gas  in  a 
plane  satisfies  a  system  of  equations  in  the  hodograph  plane  (refer¬ 
ences  1  and  2) : 


q  =  -  £°(i  -  M^)  — 

dq  p  be 

> 

dcp  _  Po  cM 
be  p  q  bq 


(1) 


Here  cp(q,0)  and  iKq,0)  are,  respectively,  the  velocity  potential 
and  stream  function;  q  and  0,  the  magnitude  and  inclination  of  the 
velocity  vector,  respectively;  p  and  M,  the  density  of  the  fluid  and 
the  local  Mach  number,  respectively;  and  pQ,  the  value  of  p  at  q  =  0. 
(See  appendix  A  for  definitions  of  all  symbols.) 

This  fundamental  system,  as  shown  previously  (reference  2),  can  be 
solved  and,  in  the  case  of  ty(q,0),  the  particular  integrals  are 


and 

6  and  JT  (1  -  t)P  (*0 

and,  in  the  case  of  cp(q^6)>  are 


(1 

(1 

(1  _  T)“0  _  iyT(l  -  t)~P  and  6  (7) 

All  these  results  are  well-known  except  equations  (3)  and  (6)  which  can 
he  found  in  appendix  B.  Here  Fv(t)  and  F_v(t)  stand,  respectively, 

for  the  hypergeometric  functions  F(av,hv;  v  +  1;  t)  and 

F^a^  -  v^v  -  v;  1  -  v;  t)  for  nonintegral  parameters  V.  In  case  v 
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is  an  integer  n,  the  definition  of  the  second  integral  is  modified. 
The  new  function  will  include  a  logarithmic  term  (cf.  reference  2)  and 
is  denoted  hy  F_n(T).  The  functions  |  (T)  and  g  (t)  are  defined 
by  . 


VMt)  =  2t  -t  lo§e  tv/2fv(t) 

V 

V|-v(t)  =  2t  ^  l0ge  t_V/2f_v(t) 


(8) 


and  the  constants  av  and  b^  are  defined  by 

av  +  bv  =  v  -  P 

aVbV  =  -  \  Pv(v  +  1) 
and  the  variable  r  is 

T  .  iaL 

2P  cQ2 

where  7  is  the  ratio  of  specific  heats  of  the  gas  and  c0,  the  speed 
of  sound  at  q  =  0. 

Furthermore ,  Fv(t)  is  an  analytic  function  of  both  variables  T 

7-1 

and  v.  For  fixed  t  <  - ,  it  has  simple  poles  at  v  =  -n  and, 

7  +  1 

accordingly,  possesses  an  expansion  of  the  form  (references  3  and  4): 


3  = 


(9) 


t  ^  \ 


Fv  (T ) 


f(i-)TV 


00 

-r 


v  +  m 


TmT(r)m+VFm(T) 


(11) 
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where 


f  (t  )  = 


-  T)rl2/' 


(1  -  T) 


(1  - 


T(t)  = 


y  -  ?i) 


[ri(l  -  r)l/2  +  (l  -  7i2t)i/2_ 

71  (1  -  t)1/2  ♦  (1  -  r^)1/2 


7i 


7  +  1 
=  \7  1 


>  (12) 


Cm  — 


r(am)r(m  +  1  -  bm) 


r(am  -  m)r(l  -  bm)ml(m  -  1)1 


This  expansion  is  valid  for  all  values  of  V  except  negative  integers 
When  V  =  -n  the  limiting  process  gives 


F  (t)  =  f  ( T )  t”^(  t  )  ~YL  '  n  +~m  TmTm  " 

-XI  * - =r  -n  -r  m. 

m=X 

cn[x%(T)  log.  T  ♦  (13) 


where  the  prime  indicates  that  the  term  m  -  n  is  to  be  excluded. 


From  these  particular  solutions  one  can  construct  a  solution  for 
any  boundary  in  the  hodograph  plane.  Once  the  functions  cp(q,6)  and 
iKq,0)  are  known,  the  flow  pattern  in  the  physical  plane  is  given  y 
integrating 


bx 

be 


Sy 

be 


p 


sin  e 


be) 


cos 


ilk) 


* 


V 
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and 


6x 

6y 

6q 


=  ifcos  0  |2  -  ^  sin  0  |i 

6q  P  ag> 

=  -/sin  e  ^2  +  £°  cos  6 

q\  oq  P 


(15) 


where  x  and  y  are  the  coordinates  of  a  point  in  the  physical  plane 


MAPPING  OF  AN  INCOMPRESSIBLE  FLOW 


The  potential  flow  of  an  incompressible  fluid  past  a  symmetric 
Joukowski  airfoil  without  circulation  is  determined  by  a  complex 
function  QW(z): 

W(z)  =  C  +  e  +  — --+  z  =  £  +  i  (16) 

£  +  e  £ 


where 


z  =  x  +  iy 


£  =  5  +  in 


Here  all  the  physical  quantities  are  made  dimensionless  by  dividing, 
respectively,  every  length  and  velocity  by  the  radius  of  the  generating 
circle  and  the  free-stream  velocity.  Then  e  is  a  purely  real  number. 
The  complex  velocity  w  for  such  a  flow  is 


£2(£  +  1  +  2e) 

U  +  1)(£  +  e)2 


This  equation  expresses  the  relationship  between  w  and  £.  It 
can  also  be  interpreted  as  a  transformation  function  from  the  £ -plane 
to  the  w-plane .  In  other  words,  if  a  flow  field  in  the  physical  plane 
is  given,  its  image  in  the  w-plane  is  immediately  known  by  means  of 
equation  (17).  Since  the  present  problem  is  formulated  in  the  w-plane 
instead  of  the  z-plane,  the  flow  field  in  the  z-plane  can  be  calculated 
from  the  inverse  transformation  £(w),  which  is  given  by  solving 
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„ 3  ,,  „  ,,2  e(e  +  2)w  €  w 

{3+(l  +  2e)t  -  i  -  — 


The  three  solutions  of  this  equation  are: 
y  /  \  \  .  \  1  +  26 


?l(v)  =  P(w)  +  Q(w) 


^2(w)  =  -||p(v)  +  Q(w)]  -  ^p|p(w)  -  Q(w)]  -  - 

i${v)  =  -|[p(w)  +  Q(w)J  +  ±p[p(v)  -  Q(w)J  -  - 


1  +  2e 

3 


where 


ij(1 .  e)3  +  Li 

3  1  1  -  v 


i33/26(l  +  26)3/2w  1/2  _  (1  -  e)3(l  +  0 

^/2 


(1  -  w) 


(1  +  26  )  3 


n  ,0  9e(l  +  €  +  e2) 

Q(w)  =  --  1(1  -  e)3  +  - ; — — - + 

3  1  -  w 


i33/2g(i  +  2g^3/2vl/2  ^  (l  -  e)3(l  +  g)  w 

(1  -  v)3/2  (1  +  26)3 


This  indicates  that  to  a  neighborhood  in  the  w-plane  there  corre¬ 
spond,  in  general,  three  distinct  neighborhoods’  in  the  £-plane  or  the 
z -plane  (see  fig.  l) .  The  exceptional  points  where  two  of  the  three 
solutions  become  equal  are  given  by  dw/d£  =  0.  This  yields  three 

points  {;  =  0,  £  =  --  +  and  £  =  °°  about  which  two  of  three 

2  +  e 
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1° in*  Furthermore,  it  remains  to  investigate  how  the  axis  of 
?'  v  =  0,  is  mapped  in  the  £ -plane.  By  separating  into 
al  and  imaginary  parts,  the  velocity  component  v  in  terms  of  $ 
and  T)  is: 

V  =  -  2eTlU2  +  fjS  +  n2)  [(2  +  e)g  +  1  +  2e~) 

Then  v  =  0  corresponds  to 

n  =  0 

(2  +  e)|  +  1  +  2e  =  0 


(5  ♦  i) 


2  2 
2  +  n2  =  e_ 


^iToTZHTnr^l  refXf  °f  the  W-plane  corresponds  to:  (a)  The 

point  parallel  to  Se^P  nS'  ^  f  ftraight  line  through  one  branch 
P  parallel  to  the  T]-axis,  and  (c)  a  circle.  Moreover  while  1, 

ncreases  from  zero  to  one  along  v  =  0,  ^(w)  decreases  from  the  for¬ 
ward  stagnation  point  to  negative  infinity,  £ (w )  increases  from  zero 
to  positive  infinity,  and  £3(v)  decreases  to  a  negative  number.  Evi¬ 
dently,  branches  I  and  II  join  smoothly  at  the  line  g  =  _  UiJl 


1  +  2e 

2  +  e 
.2 


branches  II  and  III  meet  along  the  circle  +  £.\2  +  ^2  .  e2 

Branch  III  then  lies  entirely  inside  the  circle  of  radius  e/2.  Hence 
the  flow  over  a  body  with  a  sharp  trailing  edge,  complicated  as  it  is  ' 

^  V*  S“e  as  that  °ver  ^  which  is  s»ric' 

eouat Ions  hi  ”  h  tv“  stapiation  points.  Finally,  by  substituting 
equations  (19)  in  equation  (l 6),  the  corresponding  branches  of  the  *om 
plex  potential  are  ^(w),  0W2(w),  and  pW3(v). 

Jftien  circulation  is  present,  having  the  magnitude  required  by  the 
condition  of  finite  trailing-edge  velocity  at  an  angle  of  attack  a, 
the  complex  potential  is  9 


,W(z)  -  (£  t  e).-ta  i  e1*  ,  21(1  +  o  sm  a  loge  (£  +  .)  (21) 
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with  z  =  £  +  i;  and  the  complex  velocity 


+  e)e 


The  inverse 
solution  of 


J(  £  +  e)e"ia  +  (1  ' 

U  +  DU  + 


transformation  function  £(v)  in  this  case  consists  in  the 


(e-ia  _  w)^3  +  [ee-ia  +  (1  +  e)eia  -  (1  +  2e)v\f  -  e(2  +  - 


The  solution  of  this  equation  can,  of  course,  he  given,  hut  the 
the  branch  points  are  £  -  0,  £  an 


(1  +  e 


)(l  -  e2ai)^2  +  e  (3  +  €  +  (1  +  e)e2a^  +  2€[c  +  (l  +  e 


e2ai  =  0 


for  small  values  of  a,  to  curves  in  the  £-plane. 

a(l  +  *[<2+  OS  +  l>2il'!3  +  aU+  £>[5|3  + 

6(1  +  e)i2  -  (3  -  25  -  e2)l  -  (1  +  e>!l'l2  +  eS<5  +  e)C<2  +  °*  * 


1  +  Tj  +  a(l  +  €)^  (£  +  +  e)  ^ 


t 

This  was  not  done  for  ,,  T+  should  he  remembered  that 

difficult,  if  not  impossible,  to  handle .  It  should  *e  ^  flow 

the  present  procedure  of  constructing  *  \ov  solution 

depends  primarily  on  the  knowledg  f  P  zer0  circulation,  as 

in  the  form  of  power  series.  Even  m  tne  case  u 
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can  be  seen  from  equations  (20),  the  series  expansion  of  the  complex 
potential  yields  sets  of  coefficients  which  are,  at  least,  double 
infinite  series  (cf.  appendix  C).  One  can  easily  see  that  the  initial 
task  alone  is  a  formidable  one,  not  to  mention  the  later  computation 
of  the  compressible -flow  solution.  This  amount  of  labor  is  particularly 
unjustifiable  in  this  problem  because,  even  if  one  starts  with  an  exact 
solution,  he  may  still  end  up  with  a  body  very  different  from  what  he 
expected  because  of  the  distortion  due  to  compressibility.  For  this 
reason  deviation  from  the  formal  procedure  is  recommended  wherever  it 
is  convenient.  In  this  report,  the  case  of  small  thickness  and  small 
angle  of  attack  is  studied  first. 

Before  proceeding  further,  it  is  mentioned  that  the  Joukowski  air¬ 
foil  is  characterized  by  having  a  cusp  at  the  trailing  edge..  Although 
the  velocity  is  finite  there,  the  acceleration  is  infinite  and,  con¬ 
sequently,  the  analyticity  of  the  mapping  f unction  breaks  down;  namely, 
dQw/dw  =  0.  In  order  to  preserve  this  property  in  the  case  of  compress¬ 
ible  flow,  the  stream  function  t(q ,0)  must  satisfy 


=  0 

dq 

*1=0 

Be 


(24) 


In  the  case  of  symmetric  flow,  the  first  is  satisfied  identically.  The 
second  serves  to  determine  the  speed  at  the  trailing  edge. 


APPLICATION  TO  THIN  SYMMETRIC  AIRFOILS 


Consider  the  case  of  a  thin  symmetric  airfoil  present  in  a  uniform 
flow  without  circulation,  namely,  a  =  0.  Since  the  thickness  ratio  of 
a  Joukowski  airfoil  is  directly  proportional  to  the  parameter  e ,  for 
thin  airfoils  e  must  be  small.  Suppose  that  the  airfoil  is  so  thin 
that  terms  higher  than  e  are  negligible.  Equation  (l8)  is  then 
simplified  to 


£2  +  (1  +  2e)£ 


(25) 
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1  -  2e 

where  e-i  =  - -,  being,  of  the  order  of  magnitude  unity. 

1  +  2e 

By  making  such  an  approximation  it  is  seen  that  the  third  solution 
disappears.  This  is  explained  by  the  fact  that  the  third  branch  lies 

1  2 

wholly  inside  a  circle  of  area  -jj-  jte  which,  according  to  the  present 
approximation,  is  zero;  therefore  ^j^)  =  °*  The  branch  points  are 
now  located  at  £=-i(l  +  2e)  and  £  =  «,  connected  by  the  branch 
line  (fig.  2)  2%  +  1  +  2e  =  0. 

The  problem  is  indeed  simple.  The  question,  however,  is:  What 
flow  is  represented  or,  in  other  words,  how  is  the  complex  potential 
modified?  By  using  the  simplified  transformation  function,  a  straight¬ 
forward  integration  gives 


QW(z)  =  i  +  €  +  1  |-2-€  loge  (1  +  J-) 


(27) 


The  flow  in  the  £-plane  is  then  that  produced  by  the  superposition  of 
a  sink  and  a  source  at,  respectively,  £  =  0  and  £  =  - 2e  in  a  uniform 
flow.  It  is  very  easy  to  see  that,  if  e  is  small,  the  vanishing  of 
the  stream  function  corresponds  closely  to  a  circle  with  radius  1  +  € 
and  center  at  £  =  -e.  Hence  for  small  values  of  e  the  curve  in  the 
z -plane  will  not  differ  much  from  a  Joukowski  airfoil.  By  substituting 
and  from  equations  (2 6)  in  equation  (27)  the  two  branches 

of  the  complex  potential  are 
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I 


* 


0W2(v) 


(1  + 


(28) 


1  +  2e 
2e 


log, 


1  - 


1  " 


1  +  2 


+  1 


w 


By  taking  e  =  0.08,  the  hodograph  ^(q,#)  =  0  from  equations  (28) 
and  that  of  a  symmetric  Joukowski  airfoil  are  compared  in  figure  3..  It 
is  seen  that,  except  in  the  neighborhood  of  the  singularity  w  =  1, 
the  agreement  is  quite  close  even  for  a  thickness  ratio  of  about 
10  percent. 

For  the  case  where  a  weak  circulation  is  present,  by  the  same 
method  of  approximation,  equation  (23)  simplifies  to 


s 


2 


+ 


fl  +  2e  + 


SSL..  0 

1  -  w 


(29) 


Here  a  is  assumed  to  be  the  same  order  of  magnitude  as  6,  and  there¬ 
fore  the  product  of  e  and  a  is  also  dropped.  The  solutions  of  this 
equation  are 
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ih 


the  corresponding  complex  potential  is 

0W(w)  =  (1  -  al)£  +  e  +  1  +  2*/-  —  l°ge  (x  +  j)  +  2al  l0®e  (£  +  2€) 

(31) 


The  body  in  the  £ -plane  is  again  a  circle,  provided  both  e  and  a 
are  small.  A  substitution  of  £]_  and  £2  from  equations  (30)  in 
equation  (31)  gives,  respectively,  QW1(w)  and  0W2(w) . 

It  should  be  noted  that,  in  the  previous  case  when  circulation  was 
present,  it  altered  radically  the  nature  of  the  singularity  of  the  flow. 
Now  by  comparing  equations  (30)  with  equations  (26)  one  can  see  that 
circulation  has  introduced  essentially  no  complication  and  hence  the 
problem  can  be  similarly  handled  as  in  the  case  where  no  circulation 
exists. 


EXPANSIONS  OF  QW(w) 


When  a  =  0.-  For  thin  airfoils,  as  seen  from  the  previous  section, 
the  domain  of  interest  has  two  branch  points  (;  =  -i  (1  +  2e)  and  £  = 
the  corresponding  points  in  the  w-plane  are  w  =  1  and  w  =  at 

each  of  which  GW(w)  becomes  singular.  As  a  result,  two  different 
expansions  are  required  to  represent  each  branch  of  QW(w).  In  the 
case  of  owl(v) >  the  origin  is  a  regular  point;  for  oW2(w)  it  is  a 
logarithmic  singularity. 

Consider  0W^(w)  for  the  case  a  =  0  first:  About  the  forward 
stagnation  point  | £|  ^  0  and  QW1(w)  can  be  simplified  in  the  form: 


z  \  u  1  -f-  2  € 

0Wi(w)  =  ^  +  e  +  — - - 

£l 


e(l  +  26) 

£l2 


(32) 


By  means  of  the  relation: 
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-I/2/  2  \-l/2 

owl(w)  can  he  reduced  to  two  elements  (l  -  w)  (l  -  wj  and 

(l  -  w)1/^  ^1  .  which  possess  the  following  expansions  in  the 

respective  domains: 

n(iv 


(1  -  »)l/2(l  -  ^)V2  =E  lsn<i)’'° 
x  '  0 


|w|  <  1  L  (33) 


where 


S„(1)  -  <-l)nC 

m=0 


n 


A(1)  -  (-1>° 


2(  n-m) 


2(n-m) 


L  (34) 


and. 


(1  -  -  c/w)-l/2  -  £  Sn‘°>ei2V  * 


uu 

H  sn<“>v-° 
1 


1  <  w  < 


(1  -  v)l/2(l  -  ei2v)l/2  .  £  lSn(»)ei2-V  + 

0 

£  lSn(o)v-” 


er 


(35) 


1 6 
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where 


.  .  -  /  i\/  1  \ 

;n(°)  =  (_l)n  p  2  2  , 

0  \  m  /\m  +  n / 

Sn(0>  -  <-«n  £  (2)(  2  \ 

0  \m/\m  +  n / 


A  simple  reduction  then  yields 


00 

>wl(v)  =  -Ao  -  Y^.  AnwI1  +  0(f2)  lwl  <  1 


where 


An  - 


♦  <i^i>![Bn+1<^  -  (1  ♦ 

iV/11]  *  *  ^J11  - 


'Sn+2(1)  +  (2  +  ^)Sn*l(i> 


0c  2q  (l)  q  (i)  ,  -  2q  (  1  Jl 

2el  sn  -  sn  +  €1  sn-l 


Here  Ao  is  zero,  as  is  easily  shown  from  equation  (27)j  in  the  above 

^  P 

expansion  it  would  have  been  of  the  order  of  e  and  hence  is  included 

in  the  symbol  o(e^).  The  velocity  potential  and  stream  function  are, 
neglecting 


oh  "  -Ao 


-£v" 


cos  n0 


q  <  1 


00 

=  P  Anqn  sin  n6 
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Separating  into  real  and  imaginary  parts,  Qcp1  and  ^  are, 
approximately, 


o<T>i  =  B  ^q-2  cos  2®  +  B_lq  ^ 

o\  =  B-2q”2  Sin  26  +  B-lq”1 
for  1  <  q  <  — i-. 

e-.2 


cos9  +  Bo+i:  ^Bvqv  -  Cvq“v)  sin  V9 

00  . 
sin  G  +  (BvqV  +  Cvq_V)  cos  vG 


>(42) 


In  the  case  of  QWg(w) ,  as  w  =  0  is  not  a  regular  point,  a  Taylor 
expansion  does  not  exist.  It  can  be  expanded,  however,  in  the  following 
form: 


oW2< 


w)  =  A_nw"n  +  Y1  +  °(€2)  M  <  1 


(43) 


where  An  is  defined  in  equation  (38)  while  A_0,  A_-|_,  and  A_2  are 


-o 


-2 


3(1  +  2e) 
4e 

(1  +  2g)2 

6 

(1  +  2e)3 

4e 


(44) 


On  the  other  hand,  in  the  region  1  <  |w|  <  the  expansions  are 


0w2(w: 


Y~  B  w-n  -  if  (b  wv  +  C  w-v)  +  o(e2) 

0  “n  V  v 


(45) 
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Here  the  constants  B_j,  BQ,  B  ,  and  Cy  are  defined  in  equa- 

1  tions  (4l).  It  is  clear  that  by  circling  once  about  w  =  1  branch  I 

goes  smoothly  over  into  branch  II,  as  is  easily  seen  from  equations  (40) 
and  (V?).  Similarly,  the  velocity  potential  and  stream  function  are, 
omitting  o(e2), 


oqp2  =  A_nq~n  cos  n0  +  Anqn  cos  n0 

2 


^2  =  Y~  A_nq'‘n  sin  n0  +  }  '  Anqn  sin  n0 
1  2 


0  <  q  <  1 


l  W) 


and  for  1  <  q  <  e 


-2 


0^2 


=  YZ  B-n^-n  cos  nS  +  (®vqV  -  Cvq-V)  sin 

r\  r\  \  v  • 


ve 


2  00  .  V 

1^2  =  H  B_nq-n  sin  n6  -  T~  (Bvqv  +  Cvq-V]  cos  ve 

1  ^  w 


>  (^7) 


When  a  ^  0.  -  When  a  is  finite  but  small  in  the  case  of  a  thin 
airfoil,  equations  (30)  show  that 


=  ^!(o)  +  a?!(l) 
C2  =  £2(°)  +  a^1) 


m 


where  and  are  the  transformation  functions  for  a  =  0 

and  and  are  correction  terms  due  to  circulation.  There- 

fore,  they  are  defined  as  follows: 


20 
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m 


(50) 


where  qW-^ °)  and  0W2^0^  are  defined  in  equations  (28)  and 
and  0^2  ^ are  found  to  be 


which  are  valid  if 


?(o>  1 0. 


* 
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To  expand  0 W(w)  ,  "besides  the  fundamental  forms,  expansions  (33) 

and  (35),  it  is  required  also  to  consider 


logp  (1  -  2e)(l  -  w)  +  (1  +  2e) 


(1  -  ^)1/2] 


t1  - " 


>i/2 ♦ (1  -  *a)i/2] 


For  the  first  expression  let 


F(w)  =  —  log 


dv  °e 


[(1  -  2e)(l  -  v 


1/2 

)  +  (1  +  2e 


)(i  -  e/w)1^ 


By  direct  differentiation  one  obtains 


F(w)  =  (1  -  w)  1/2(±  - 


By  making  use  of  equations  (33)  and  (35)>  a  term-by-term  integration 
gives 


loge  |(1  -  2e)(l  -  w)l/2  +  (1  +  26 ) (l  -  6 


el  t£2-sn_n(i) 

- r -  v11  +  l°Se  2 

2  ±  n  Be 

and,  aside  from  a  constant, 


M  <  1  (52a) 


[u  -  2e)(l  -  v)1/2  +  (1  +  2e)(l  -  e^v)1^  = 

^  H-Sn+1(0)-“V  +  Sn(°)612nwv) 


(52b) 
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A  similar  method  yields 


[(1  -  v)1/2  +  (l  -  e^w)1^  = 


°2L  G„  \  —  / 

-  i  7^  — -  v11  +  log  2 

p  £-=—  n  e 


M  <  1 


(53a) 


>ge  [(1  -  v)1/2  +  (l  -  e-^v)1/2]  = 

|k  i(sn+i(o)«i2,1+2»v  -  Sn<0)»-1)  1  <  1*1  <  «l'8  <53l>) 


Furthermore  y 


♦  0 


can  be  written  in  an  expansible  form: 


“  — ll  +  € 


II  -  w 
1  -  e^w 


Collecting  all  the  individual  expansions,  one  finds 


-i  =  2«i  + 


+  o(e2) 


w|  <  1 


where 


A„(l)  ‘  -(l  ♦  SnU))  *  ^-^(Sn(1)  -  el2sn-lU)) 

HriH(1-‘i2s"-i(i))-^»(1-s“U)) 


1-26 

2e 


(sn-lU)  -  sn(1)) 


(55) 
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The  corresponding  functions  and  are 


®  (D  =  (*  - 0) 

o\L  he 


-  £  (-  - 
^Bv^l)qv  -  Cv^1^q"v)  sin  v6  + 


sin  n0  + 


Constant 


oh 


(1) 


00  . 

=  -1  *■  h  loge  q  +  5Z  +  ^)q_n  COS  n0  + 

^  (Bv(l)qv  +  Cv(l)q_V)  cos  V0  +  Constant 


Similarly,  for  the  second  branch  the  expansion  is" 


-i  W  =  C  rti 

o  2  o  2e  2e 


—  loge  w  +  YL  Cjj^h11  +  o(e2) 


where 


(l)  1  +  he  ,  1  n  1  +  2e 

CD^  '  =  — - -  loge  e - log  — - -  +  2 

0  he  e  2e  e  2 


cn(l)  =  -(1  -  S^)  -  ^p(sn(i)  - 

(l  +  ^lsn-l(l))  ~  ^  +  Sn(l))  + 


1+4 

4en 


1  +  2e 
2€ 


(sn(l)  -  Sn-l(i))6l 


>  (59) 


M  <  1 

(60) 


> 


(61) 
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-2  . 


For  the  region  1  <  |w|  <  ,  it  is 


00  .  » 

-i  ow2(l)  =  lo8e  w  -  i°ge  veJtl  +  XZ  ^  +  5yw-n  “ 

y2"  ^Bv^1^wv  +  cvw"v)  +  o(e2) 


(62) 


where  a  constant  of  order  unity  is  again  left  out.  The  corresponding 
real  and  imaginary  parts  are: 


oCp2(l)  =  ^  (*  -  0)  +  IZ  Cn(l)qn  sin 


n0 


q<!  >  (63) 


00 

^1'1  =  Co(l)  ■  f**q+r  Cn^^qn  cos  n0 

hi  6  O 


-2 


and,  aside  from  a  constant ,  in  1  <  q  <  , 


=  ^ +  ^1  +  sin  ~ 

II  (Bv(l)qV  -  Cv(l)q“V)sin  V0 


“I 


>  (64) 
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CONSTRUCTION  OF  A  SOLUTION  FOR  A  SYMMETRIC  AIRFOIL 

The  complex  potential  0W  for  a  flow  over  a  tllin  JQui;owskl  airfoil 
has  now  heen  transformed  and  expanded  in  power  series  of  w  for  both 
the  case  of  zero  circulation  and  of  a  weak  circulation.  Knowing  the 
form  of  expansion  of  qW>  a  solution  relating  to  a  flow  of  incompres¬ 
sible  fluid  about  a  certain  body  can  be  obtained  immediately  by 
replacing  each  term  of  the  series  by  the  proper  particular  integral  as 
given  in  the  section  "Fundamental  Equations  and  Their  Particular 
Solutions."  For  instance,  the  complex  potential  in  the  case  of  zero 
circulation® is 

00 

DWi(w)  =  -  Anwn  |v|  <  1 
0 

The  complex  potential  Wp(w/r)  relating  to  a  compressible  flow  is, 
accordingly, 

Wi(»,T)  -  -  Y1  V„(r)(T)v'>  ♦  ti  Vm(fiT)"  ' 

0  m=2 

w<i  (65) 

of  which  the  imaginary  part  gives  the  stream  function,  namely, 

tyl(q,e)  =  lnijw1(w,-0j  (66) 

Here  the  definition  of  Fn^r^(T)  may  be  either  Fn( tJ/f^Tl)  (refer¬ 
ences  1  and  2)  or  Fn(r)/T1(r1)  (references  3  and  4),  where  t  *  t± 
when  q  =  1.  It  follows  from  equation  (66)  that  the  stream  function  is 

ii_(q,®)  =  Y~  AI1qnFn^ r ^ ( t )  sin  n6  +  ^  hmcm(friTl)  V  sin  m6  ~ 

2  2 

^2  sln  nl-2 


q  <  1  ( 67 ) 
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From  equation  (5)  the  velocity  potential  is 

00 

=  Aq  -  (l  -  t)~P  AnqnFn( r  ^  ( T )  ln( T )  cos  ne  + 

2 

^  hmcm(friTl)mV^m  cos  -  A2c2  ^  [(fr1T1)n1|(n5n  gos  n©j  ^ 

q  <  1  (63) 

Both  series  (equations  (67)  and  (68))  converge  only  within  the  unit 
circle  q  =  1.  In  order  to  have  a  function  to  cover  the  whole  field  of 
flow*  both  series  must  be  continued  across  the  circle  of  convergence. 
There  have  been  two  alternative  methods  proposed  to  effect  this  con¬ 
tinuation.  The  first  method  (references  1  and  2)  is  to  have  the  solu¬ 
tion  corresponding  to  the  flow  in  the  upper  z-plane  arranged  in  one 
w-plane.  After  the  series,  for  example,  equations  (39),  (42),  and  (46), 
are  modified,  they  are  required  to  he  continuous  with  continuous  normal 
derivatives  on  the  circle  q  =  1  (fig.  4).  On  account  of  compressi¬ 
bility,  the  modified  solution  will  he  discontinuous  in  normal  derivative 
and  hence  W-j__  and  W-j_  can  no  longer  represent  the  same  function.  To 

correct  this,  one  set  of  coefficients  is  left  free  and  an  extra  series 
is  added  to  W1*  These  two  sets  of  coefficients  can  then  be  determined 

hy  the  condition  at  q  =  1  (reference  2).  The  second  procedure  (refer¬ 
ences  3  and  4)  is  to  proceed  from  a  Taylor  expansion,  equation  (67),  say. 
Expressing  the  function  Fn(-r)  in  the  form  of  equation  (ll)  gives  rise 
to  a  convergent  double  series  which,  when  the  order  of  summation  is 
changed,  is  transformed  into  a  form  being  convergent  even  for  values  of  q 
greater  than  unity.  Analytical  continuation  is  thus  accomplished.  In  the 
following  discussion,  the  latter  procedure  is  adopted  for  simplicity. 

At  this  point,  the  convention  for  numbering  the  branch  of  0W(w) 
is  altered. .  Instead  of  dividing  the  z-plane  fore  and  aft,  it  will  be 
more  convenient  to  divide  the  plane  into  upper  and  lower  halves.  Then 
owl(w)  indicates  the  upper  half  and  QW2(w),  the  lower.  The  parts 

corresponding  to  the  regions  q  <  1  fore,  1  <  q  <  €l"2,  and  0  <  q  <  1 
aft  are,  respectively,  QV1_)  Qvlt  and  QW1+. 


By  the  second  method,  replacing  Fn(r)  in  equation  (65)  by  equa¬ 
tion  (11)  and,  for  evident  reasons,  adopting  the  definition 
(r)  / 

Fn  ( T)  -  Fn(TVTl(Tl)>  equation  (65)  can  be  written  as 


=  -f(T)  y  An(tw) 


An(tw) 


TmTmFm( 
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where  t  =  t/Tj_.  The  first  series  can  he  summed  and,  hy  changing  the 
order  of  summation,  the  second  series  can  be  transformed  into  a  single 
series.  Thus 

00  Pw 

W!(w,t)  =  f(r)  ^(w)  -  Yi  Cm(^F1T1)m1l;m(T)eim0  "“"Vi  dw  + 

m=l  '  UO 

£  ^  <© 


where  w  =  tw  and  ^(t)  =  r  Fm(x).  The  first  term  corresponds  to 

that  resulting  from  the  asymptotic  summation  (references  1  and  2)  and, 
in  the  special  case  7  =  — 1,  reduces  to  the  Karman— Tsien  approximation, 
The  integral  in  the  second  term  is  bounded  and  the  series  can  be  shown 
to  be  convergent  for  q  >  1  (references  3  and  4). 

By  taking  a  path  starting  from  the  forward  stagnation  point  to  an 
arbitrary  point  in  the  annular  region,  namely, 


w  PV 1  r« 

wm“d0Wj_( w)  dw  =  I  wm_dQW1(w)  dw  +  I  w111-  0W-]_(w)  dw 


where  3^  is  a  fixed  point  in  the  annular  region,  a  simple  transforma¬ 


tion  leads  to 


.  fz  B_nF_n(r>(T)„-“  +  i  fz  *  <y?.v<r>(T)v-v] 


■k2c2  +  b.2c2  H- 
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by  replacing  QW1  by  the  proper  expansion  under  the  integral  sign, 
where 


hm  - 


m-1 


ow! dw  -  n 


'V 

B-nwl 


V  +  m 


^rw| 

-V  +  m  / 


B_2  loge 


(71) 


The  integral  in  equation  (71)  can  be  evaluated  in  the  £ -plane;  namely. 


wi 

*"-Vi 

Jo 


dw 


w"'Vi  f  at 


The  fact  that 

on 

Equation  (70),  evidently,  is  the  continuation  of  equation  (65)*  Pro¬ 
ceeding  further  to  a  fixed  point  Wg  inside  the  unit  circle  hut  below 

the  real  axis,  by  describing  a  loop  surrounding  the  point  w  =  1,  the 
following  proper  transformation  is  arrived  at: 


"(m)nvine] 


is  a  solution  is  shown  in  appendix  B. 


W1  =  II  Vm^irv^  +  n  A_2F_n(r)(T)  v-n  Hr 


^  AnFn(r)( 


t)^ 


(72) 


where 


^(£2)  ■ 


2  .  ~  -n+m  „  .  -v  n+m 


^o 


2  v“-l QWi  ^  d£  - 

°  1  d£  -n  +  m  *q~  n  +  m 


A-2  loge  w2 


(73) 
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The  path  is  quite  arbitrary;  in  order  to  preserve  the  property  of  sym¬ 
metry  of  the  function  it  is  so  chosen  that  the  hj^'  s  are  real.  Equa¬ 
tions  (65),  (TO),  and  (72)  together  represent  the  first  branch  of  the 
Riemann  surface  covering  the  entire  upper  z-plane.  Separating  into 
real  and  imaginary  parts,  the  stream  function  is 


00  - 1 

tl(q,0)  =  ^|BvqvFv(r)(T)  +  Cvq_vF_v( r ) ( t)J  cos  V0 


1  <  q  <  €■ 


-2 


(7^) 


and 


00 

*l+(<l,e)  =  YZ  v(^2)cm(fiT])m^T)  sinme  - 

2 


AnqnFn^r)  ( t)  sin  n0  + 

2 


sin  n0 


<1  <  1 


(75) 
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The  velocity  potential  corresponding  to  equations  (7^-)  and  (75)  is 


cppUjG)  =  (1  -  t)"P  J 


YL  B-n^'nF-nr)5-n  cos 


nQ 


fl  +  Cvq-vF_v(r)(T)|_J  sinve 


(A2c2  '  B-2c2)^[(J7lTl)n^n^n  cos  n®jn=. 


1  <  <1  <  «i 


-2 


and 


9i+(q,0)  =  (1  -  t)“P«< 


-  ^  hffl^2)Cm  (]f7lIFl)m'l/ni' T )  cos  - 


Y~  A_nq"nF_n(r)(-r)  |_n  cos  n0  + 


Fn(r)(T)ln 


cos  n0 


(76) 


0  <  q  <  1  (77) 
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By  taking  a  path  from  the  forward  stagnation  point  hut  in  the 
reverse  direction,  the  second  branch  of  W(w,t)  can  he  joined*  Nov 
the  solution  QW(v)  ^ 0  considered  in  the  integral  appearing  in 
equation  (69)  will  he 


oW2-  =  ‘  E 
0 

0W2  =  Y1  B-nw_n  - 1  E  (vv  +  Cvw"V) 

0  O' 

2  00 

oW2+  =  E  A-nw"n  +  E  V+n 

1  0 


M  < 1 


1  <  |w|  <  e^2 


(78) 


tw|  <C  i  *  j 


Complete  symmetry  with  respect  to  the  u-axis  makes  it  unnecessary  to 
repeat  the  process. 

In  the  case  of  the  flow  with  circulation,  the  complex  potential 
consists  of  two  parts  as  shown  in  equations  (50).  The  one  which  corre¬ 
sponds  to  zero  circulation  is  given  in  equations  (37)>  (4o),  and  (43). 

The  functions  for  the  compressible  flow  are  then  given  from  equations  (67) 
to  (77).  The  part  which  takes  into  account  the  effect  of  a  weak  circu¬ 
lation  is 


1  owl-(l)  =  E  An(l)vn 
0 


|w|  <  1 


-i  .WiU)  -  B0  -  ±  **  ♦  £  (l  ♦£)^>  +  £  (bvU) 


vv  + 


cv<Dv-v)  - 


>  (79) 


loge  ve 


Ttl 


00 

- 1  oWl+(l)  =  lo®e  w  +  E  cn(l)wI1  0  <  |w|  <  1 
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By  the  convention  adopted,  qW^1),  ,  and  QW1+(l)  form  the 

first  "branch  of  and  correspond  to  the  whole  upper  z-plane.  The 

same  method  yields  a  solution  for  a  compressible  flow: 
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In  a  similar  fashion,  the  second  branch  can  also  be  connected  by- 
taking  a  path  in  the  reverse  direction  and  introducing  successively  the 
functions: 


i  ow2-(l)  =  fl  An(l)wI1 


|w|  <  1 


i  ow2 


(1)  1  +  4e 


t-ii  log,  we"1  +  (l  +  ;)»'n  - 


£( 


Bv^-*-^wv  +  Cv^^w”v)  +  Constant  1  <  |w|  < 


-i  QW2+(1)  =  -i-  logg  w  +  Cn(l)wn  |v|<l 

n 


(83) 


However,  it  should  be  added  that,  since  the  functions  W1^1^  and 

do  not  merely  differ  by  a  sign,  the  hm^^'s  in  W2+^^  and 

will  not  generally  be  equal.  As  a  result,  the  two  branches  may  not  be 
joined  with  each  other  across  the  branch  line  y  =  0,  ^  <  x  <  00.  In 

order  to  avoid  this  difficulty,  the  lower  limit  in  the  integral,  for 
instance,  equation  (73) ,  should  be  adjusted  to  make  them  equal. 


Having  the  stream  function  and  the  velocity  potential  so  chosen, 
the  flow  pattern  in  the  physical  plane  can  be  calculated  by  integrating 
equations  (l^)  and  (15)  to  give  the  coordinate  functions  x(q,0)  and 
y(q,6).  When  the  flow  is  free  of  circulation,  the  stream  function  is 
antisymmetrical  about  ‘the  x-axis.  By  taking  x  =  x-^,  y  =  0  as  the 
stagnation  point  and  x  =  y  =  0  as  the  trailing  edge,  the  constants 

of  integration  can  be  determined  on  the  unit  circle  for  6  >  0.  By 
symmetry,  the  reversal  of  the  sign  of  6  gives  the  lower  half  of  the 
flow. 


When  circulation  is  present,  the  symmetry  property  is  destroyed 
and  hence  both  branches  have  to  be  calculated  separately.  Since 
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both  w(°)  and  are  properly  joined  functions,  the  determination 

of  integration  constants  on  the  unit  circle  for  both  branches  vould  be 
sure  to  give  a  pair  of  continuous  coordinate  functions. 


Cornell  University 

Ithaca,  N.  Y.,  August  10,  19^3 
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Fv(t),F_v(t) 


t  =  i.  £_ 

2p  Co2 


7-1 


!v(t),S_v(t) 
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APPENDIX  A 

SYMBOLS 

magnitude  of  velocity  vector 
inclination  of  velocity  vector 
value  of  p  at  q  =  0 
density  of  fluid 
local  Mach  number 
stream  function 
velocity  potential 
hypergeometric  functions 
nonintegral  parameter 
integers 

parameters  of  hypergeometric  functions 

functions  defined  by  equations  (8) 

ratio  of  specific  heats  of  gas 
speed  of  sound  at  q  =  0 
defined  by  equations  (12) 
defined  by  equations  (12) 
defined  by  equations  (12) 
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x,y 

0w(z) 

z  =  X  +  iy 
e 

v 

P(w),Q(w) 

u,v 

£(v) 

Ci(w),£2(w),53(w) 

0Wi(w),oW2(w),ow3(w) 

a 

c  _  1  -  2e 

t, - 

1  +  2€ 

0W(W) 

*0 

tPo 

Sn(o),iSn(o) 

An 

o'Pl'O^l 


coordinates  of  point  in  physical  plane 
complex  potential 

geometric  parameter  of  "body 

complex  velocity  defined  hy  equation  (17) 

defined  by  equations  (20) 

velocity  components 

transformation  function  defined  by 
equation  (l8) 

solutions  of  equation  (l8)  given  by  equa¬ 
tions  (19) 

branches  of  complex  potential 
angle  of  attack 


complex  potential  for  incompressible  flow 
in  w 

stream  function  for  incompressible  flow 
velocity  potential  for  incompressible  flow 

defined  by  equations  (3*0 

defined  by  equations  (36) 

coefficient  defined  by  equation  (38) 
defined  by  equations  (39) 


4o 
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0^2 >0^2 

t  (o)  f  (o) 

*1  >  *2 

t  (1)  t  (1) 

S>1  f  =2 

Tj  (  O  )  y  (o) 

owl  >o"2 

y  (l)  yc( l) 
o"i  >o"2 


B  ^  C  A) 
•»V  >UV 


defined  By  equations  (4l) 
defined  "by  equations  (46)  and  (47) 

transformation  functions 

correction  terms  due  to  circulation 
defined  By  equations  (28) 

defined  By  equations  (5l) 
defined  By  equations  (58) 


0ep1(l),0i|r1(l)  defined  By  equations  (59) 

co^1^cn^1^  defined  By  equations  (6l) 


o^l)  >0*2(1)  defined  By  equations  (63)  and  (64) 

W1(w,t)  defined  By  equation  (65) 


Fn(r)(r)  = 


T1 

T1 

*1- >*1>*1+ 

CPl.^CPl^CPl+ 


value  of  t  at  q  *  1 
value  of  T  at  q  =  1 
tranches  of  stream  function 
"branches  of  velocity  potential 
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owl(*) 

defined  by  equation  (69) 

> 

-p 

II 

^mC7)  =  *r1B/2Fm(T) 

hm 

defined  by  equation  (7l) 

o‘,l-(l).o>'l(l).o»l+<1) 

correction  terms  due  to  circulation 

defined  by  equations  (8l) 

defined  by  equations  (82) 

o',2-<1)>o»2(1)»o«2+<1) 

defined  by  equations  (83) 
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APPENDIX  B 

PARTICULAR  SOLUTIONS  OF  FUNDAMENTAL  EQUATIONS  OF  MOTION 

jf  -'-n®  is  a  particular  integral  of  differential  equations  ( 1 ) , 

-ftne1”6)  is  also  a  solution.  For  the  differential  equation  then 


becomes 


(?  >■)*.  ■ » 


As  i|rn(q)  satisfies  the  first  equation,  the  second  equation  determines 
di| rnj  bn.  Now 


Pq 


dq)  =  J£(l  -  M2)  ||  dq  +  ^  q 


P  ’  bq 


By  substituting  ^fn(<l)eine3  for  is  eas^  to  show 

lb/.  po<*  ^n  +  < 


X  U  I  •  u  _ Ai.  1  w  _|_  _  .  - 

^  n  Sql  1  p  3nSq  p  dq  n  P  Sq 


1  b_(_±  W  +  e  £05  +  i  £°i  5nV in6 

n  5e  r1  P  dndq  P  Sq  n  p  bq  / 


Therefore 


cp  =  -1 


'o l  b_l 1  eine) 

p  5n\n  dq  ✓ 
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APPENDIX  C 


SERIES  EXPANSION  OF  COMPLEX  POTENTIAL 


The  Taylor  expansion  of  the  complex  potential  of  a  flow  over  a 
symmetric  Joukowski  airfoil  is,  for  0  <  |w|  <1, 


OO 

owl  =  An”n 

oW2=-!t^-f£B^v 


0W3  “  2 


where 


An  =  Cn  +  (1  +  «)  Dn 


Bu  =  Cv  +  (l  +  «)  D. 


gu  1  £ll  f-  f-  AY J  YkV  -J  -  *  V-i)n'332J'k  —  *  c  -  e2)J~gkU  - 


k=0  1=0  \j 


2k/W\n  +  k  -  j  -  1} 


^  1 1 1  ®U  i  Xk ;  9C  - : ; :  l  ■* 


(1  ♦  2E)3(J-k+I) 

i  V  .'P-J.2J-K-4  €J(l  +  c  +  e2)J~2ktl(l  -  <03; (1  +  c)1 


+  c)J+2(l  -  e)3V  +  2<)3(k’2)(l  +  C  +  e2)J' 


2(1  +  2e) 


--  U  +  g)3+Z(l  -  g)3Z(l  +  2c)3^k  l+2\l  +  €  +  e2)^  ' 


V  =  n  +  | 
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